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(57)Abstract: 

PROBLEM TO BE SOLVED: To provide graphical modeling method and device, 
capable of quickly and precisely inferring a gene network from a manifestation profile. 
SOLUTION: The method includes steps in which a gene manifestation profile is 
inputted, correlation coefficient matrix for genes is calculated, cluster analysis is 
executed, a cluster analysis result is outputted. and the correlation coefficient matrix 
of representative genes is calculated. Then transformation into partial correlation 
coefficient matrix is executed, minimum elements are retrieved, inverse 
transformation into the correlation coefficient matrix is executed, and significance is 
determined. When there is no significance, as a result of significance judgment, 
operation is returned to the transformation into the partial correlation coefficient 
matrix, but when there is significance, a partial correlation coefficient matrix is 
outputted, to execute graphical modeling for inferring the gene network. 
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Japanese Laid-Open Patent Publication 
N . 175305/2002 {Tokukai 2002-175305) 

A. Relevance of the Above-identified Document 

The following is a partial English translation of exemplary 
portions of non-English language information that may be 
relevant to the issue of patentability of the claims of the present 
application, 

B. Translation of the Relevant Passages of the Document 

See also the attached English Abstract. 

[INDUSTRIAL FIELD OF THE INVENTION] 

The present invention relates to a graphical modeling 
method for estimating a gene network, and an apparatus used 
therefor. 

[0002] 
[PRIOR ART] 

With the recent advance in the DNA micro array technique, 
it is now possible to simultaneously measure expression of 
thousands of genes under different conditions, which include 
different stages of cell cycle, different patterns of cell 
differentiation, diversity of somatic cells, different clinical 
conditions, and different species. Some of these different 
conditions can be ordered. For example, different stages of cell 
cycle or cell differentiation can be ordered on a time scale. 
[0003] 

As used herein, the "gene expression profile'' is used to 
refer to the expression level of a group of genes measured under 
different conditions. By evaluating patterns of expression profile, 
important information can be obtained concerning gene 
functions and regulation mechanisms. For example, there has 
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been developed a method of cluster analysis, in which genes on 
a microarray are clustered based on similarities of expression 
profiles. As used herein, the "cluster analysis'' refers to the 
process in which groups of genes that show similar expression 
patterns on a microarray under different measurement 
conditions are identified and sorted in clusters. For example, 
the hierarchical clustering (Eisen et al. 1998), the 
self-organizing mapping (Tomaya et al. 1999), and other 
clustering methods (Ben-Dor et al. 1999) have been applied to 
expression profile data. 
[0004] 

Gene clustering in expression profiling is useful in 
predicting functions of gene products, or identifying groups of 
genes that are regulated by the same mechanism. 
[0005] 

One of the important information contained in the gene 
expression profile is the expression network of genes. The 
expression level of a gene is controlled by other genes, either 
directly or indirectly. Here, such a network of genes is referred 
to as a "gene network. 
[0006] 

Estimating a gene network from expression profiles is one 
of the important goals of functional genomics. The following 
methods have been addressed for the estimation of a gene 
network: Boolean network (Somogyi and Shiegoshi, 1996), a 
differential equation method (Chen et al., 1999; D'haeseleer et 
aL, 1999), and modeling using a combination of these methods 
(Akutsu et al., 1998). 

[0007] 

[PROBLEMS TO BE SOLVED BY THE INVENTION] 
Conventional methods include estimation using 
differential equations in combination with gene disruption 
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experiments, and modeling based on Boolean modeling. 
However, neither of these methods has been put to actual 
application, 

[0008] 

The present invention was made in view of the foregoing 
problems, and an object of the present invention is to provide a 
graphical modeling method that can quickly and accurately 
estimate a gene network from expression profiles. The invention 
also provides an apparatus for such a graphical modeling 
method. 

[0009] 

[MEANS TO SOLVE THE PROBLEMS] 

In order to achieve the foregoing objects, the present 
invention provides (1) a graphical modeling method for 
estimating a gene network, the method including the steps of: 
entering a gene expression profile; calculating a gene correlation 
coefficient matrix; performing a cluster analysis; outputting a 
result of cluster analysis; performing calculations for a 
correlation coefficient matrix of representative genes; converting 
the correlation coefficient matrix into a partial correlation 
coefficient matrix; searching a minimum element; re-converting 
the partial correlation coefficient matrix into a correlation 
coefficient matrix; determining a significance; and returning to 
step (f) if there is no significance, and outputting the partial 
correlation coefficient matrix if there is a significance. 
[0010] 

In order to achieve the foregoing objects, the present 
invention provides (2) a graphical modeling apparatus for 
estimating a gene network, the apparatus including: an input 
section for entering a gene expression profile; a calculation 
section for calculating a gene correlation coefficient matrix; a 
cluster analyzing section; a calculation section for calculating a 
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correlation coefficient matrix of representative genes; a 
converting section for converting the correlation coefficient 
matrix into a partial correlation coefficient matrix; a searching 
section for searching a minimum element; a re-converting 
section for converting the partial correlation coefficient matrix 
into a correlation coefficient matrix; a determining section for 
determining a significance; and an output section for outputting 
the partial correlation coefficient matrix, the converting section 
being operated if there is no significance, and the output 
section outputting the partial correlation coefficient matrix if 
there is a significance. 
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